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ARSTRACT 


The  major  elements  of  an  approach  to  the  development  of  a  more  powerful 
spray  combustion  model  is  discussed.  A  model  is  sought  that  treats,  in  a 
numerically-ef f i 'lent  manner,  unsteady,  multidimensional,  turbulent  reacting 
flows  with  polydisperse,  multicomponent  sprays  and  properly  accounts  for 
transient  droplet  heating  and  vaporization,  liquid-phase  mass  diffusion,  and 
droplet  draz  and  trajectory.  Major  advances  are  reviewed  and  important  issues 
remaining  for  resolution  are  discussed. 


i\TTR^nrcTtn\- 


There  has  been  substantial  progress  in  the  theory  of  fuel  dronlet 
vaporization  and  in  sprav  combustion  theory  over  the  past  decade.  Reviews  of 
the  research  can  be  found  in  References  1  through  5. 

It  is  appropriate  now  to  begin  to  develop  computer  codes  that  can  predict 
the  behaviors  of  spray  combustion  flows  where  mutlicomponent ,  polvdisperse 
liauid  fuel  sprays  are  burned.  In  particular,  steady  and  unsteady 
multidimensional,  recirculating  flows  must  be  analyzed.  Account  must  be  taken 
of  the  details  of  the  transient  heating  of  individual  fuel  droplets  and  of  the 
fluid  dynamics  of  the  droplet  interior  and  of  the  surrounding  gas  film  since 
it  affects  the  droplet  heating  and  vaporization. 

Resolution  on  a  scale  of  a  few  microns  may  be  required  in  order  to 
resolve  droplet  behavior.  Olobal  dimensions  of  the  combustor  are  typically  of 
the  order  of  ten  to  one  hundred  centimeters.  This  implies  that  the  disparity 
of  length  scales  can  be  five  to  six  orders  of  magnitude.  It  is  not  possible 
therefore  to  resolve  the  problem  by  exact  numerical  solution.  Rather,  the 
droplet  behavior  must  be  modeled;  essentially  a  sub-grid  model  is  required. 
Even  though  only  representative  droplets  in  any  neighborhood  are  to  be 
considered,  hundreds  or  even  a  thousand  of  these  average  droplets  must  be 
analyzed.  Computational  simplicity  and  speed  are  extremely  important  issues 
which  must  be  addressed. 

It  is  the  purpose  of  this  paper  to  indicate  and  to  critinue  some  of  the 
kev  elements  in  the  development  of  a  physical lv-accurate  and  numerical lv- 
efficient  computer  code.  Major  advancements  will  be  reviewed  and  remaininp 
issues  for  further  resolution  will  be  outlined.  Other  spray  models  have,  of 
course,  been  developed  in  the  past  and  are  well  reviewed  ^  The 


discussion  herein  Is  significantly  different  in  that  an  approach  is  needed 
that  treats,  in  a  numericallv  efficient  manner,  unsteady,  multidimensional 
(including  recirculating),  turbulent  reacting  flows  with  multicomponent, 
polvdisperse  spravs  and  properly  accounts  for  transient  droplet  heating  and 
vaporization,  liquid-phase  mass-diffusion,  and  droplet  drag  and  trajectory. 

Ir  the  next  section,  a  simplified  droplet  model  is  discussed.  In  order 
to  calibrate  that  simplified  model,  an  exact  solution  is  required.  Also, 
droplet  drag  coefficients  must  be  determined  before  good  trajectory 
predictions  will  result.  These  issues  are  discussed  in  the  third  section. 
Sprav  combustion  in  one-dimensional,  two-dimensional,  and  axivsrmetric 
configurations  are  examined  in  the  fourth  section.  Some  concluding  comments 
are  given  in  the  final  section. 

SIMPLIFIED  ALGORITHM  FOR  PREDICTION  OF  DROPLET  VAPORIZATION' 

The  fuel  droplets  in  a  combustion  chamber  are  formed  during  tbe 
atomization  process  and  move  through  an  environment  of  hot  gases.  The 
magnitude  of  the  relative  velocity  between  the  droplet  and  the  gas  implies 
that  the  Revnolds  number,  at  least  initially,  is  laree  compared  to  unity.  A 
laminar  boundarv  layer  exists  at  the  surface  of  the  droplet  in  the  gas  moving 
past  the  droplet.  The  shear  force  at  the  droplet  surface  causes  an  internal 
circulation  of  the  liquid  in  the  droplet,  typically  with  a  boundary  layer  at 
the  surface,  an  internal  wake  along  the  symmetry  axis  and  a  toroidal,  nearlv 
inviscid  vortical  core  (see  Figure  1).  Tbe  fluid  flow  is  established  verv 
quicklv  and  can  be  considered  as  quasi-steadv. 

The  high  temperature  environment  means  that  the  droplet  vaporizes  quicklv 
and  the  droplet  lifetime  is  comparable  to  the  time  required  for  a  thermal  wave 


to  penetrate  to  the  droplet  Interior.  This  rapid  regression  results  in  the 
highly  transient  nature  of  the  problem.  The  internal  circulation  causes 
significant  convective  heating  of  the  interior;  conduction  is  still  dominant 
in  the  direction  normal  to  the  liquid  st reamsurf aces . 

The  earlier  theories^ * 2 ^  for  convective  transient  heating  and 
vaporization  involved  treatment  of  the  auasi-steadv  laminar  axisvmmetric 
boundary  layers  in  the  gas  and  in  the  liquid  at  the  surface  by  solving  the 
governing  two-dimensional  partial  differential  equations  by  integral  methods 
and  numerical  analysis.  The  vortex  core  analysis  is  simplified  by  the 
approximation  that  streamsurf aces  are  surfaces  of  uniform  temperature  and 
composition,  although  the  values  are  time-varying.  This  high  Peclet  number 
approximation  allows  heat  and  mass  diffusion  only  normal  to  streamsurf aces .  A 
coordinate  transformation  to  streamsurface  coordinate  yields  a  reduction  of 
the  axisvmmetric,  unsteady  core  diffusion  problem  to  a  one-dimensional, 
unsteady  problem. 

These  earlier  theories  demonstrated  that  transient  heating  and 
vaporization  and  internal  circulation  are  of  critical  importance  and  showed 
many  interesting  details  about  the  phenomena.  However,  they  were  costly  in 
terms  of  computer  time.  Most  of  the  computer  time  was  related  to  the  two- 
dimensional  resolution  of  the  gas-phase  boundary  layer  and  of  the  liquid-phase 
boundary  layer  at  the  droplet  surface.  More  recent  efforts  by  Tong  and 
Sirignano  have  resulted  in  a  simplified  algorithm  for  the  calculation  of 
transient  droplet  and  vaporization.  Through  the  use  of  approximate  similar 
solutions,  the  gas-phase  boundary  layer  was  analyzed  by  a  one-dimensional 
method^9,10,1 1 ,12^.  In  another  analysis/1^  it  was  shown  that  certain 
details  of  the  liquid  boundary  layer  profile  could  be  neglected,  resulting  in 
further  simplifications. 


The  simplified  formulation  assumes  a  Blasius  profile  in  the  gas-phase 


houndarv  laver  governed  by  the  following  third  order-equation 
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The  boundary  conditions  are 


(i)  at  n  =  0,  ^  =  0 
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The  heat  flux  into  the  liquid  is  given  by 
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The  mass  flux  due  to  vaporization  is 


The  modified  transfer  number  can  be  represented  by 
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Also,  the  velocity  component  in  the  major  flow  direction  is  proportional  to 
the  first  derivative  of  the  Blasius  function  f.  n  is  the  similarity 


variable,  A  is  a  constant  weighted  over  the  spherical  surface,  h  is  enthalpv, 
T  is  temperature,  k  is  thermal  conductivity,  is  fuel  vapor  mass  fraction,  v 
is  the  velocity  component  normal  to  the  surface,  o  is  density,  and  L  is  latent 
heat  of  vaporization.  The  subscripts  s,  and  e  refer  to  liquid,  surface, 

and  gas  environment,  respectively. 

Tn  standard  fashion,  the  temperature  and  concentration  profiles  through 
the  gas-phase  boundary  layer  can  be  related  to  the  velocity  profile  or,  in 
other  words,  to  the  derivative  of  the  Blasius  function.  So  solution  of  the 
Blasius  equation  will  quicklv  lead  to  solutions  of  the  scalar  properties. 
Note  that  a  phase  equilibrium  condition  relates  Y^g  directly  to  Tg . 

By  numerical  integration  of  the  above  equation,  the  Blasius  function 

f(n,B’)  is  determined  at  the  liquid-gas  interface  as  a  function  of  the 

modified  transfer  number  B'.  The  mass  vaporization  rate  per  unit  surface  area 
is  directlv  proportional  to  f(0,B'). 

The  particular  value  of  B'  depends  upon  heat  flux  into  the  liquid 

interior  and  surface  temperature  so  that  we  see  a  strong  coupling  between  the 
liquid  and  gas  phase  phenomena.  We  also  begin  to  see  the  dangers  of  many 
analyses  that  neglect  liquid  heating  or  make  assumptions  about  the  surface 
temperature.  A  useful  theory  must  allow  the  governing  physics  to  determine 
fully  the  time-dependent  value  of  B*  and  must  not  force  some  unnatural  value 
of  B'  into  the  calculation. 

The  gas-phase  analysis  has  now  produced  a  function  f(0,B')  that  can  be 
tabulated  for  a  range  of  B’  values.  Through  the  use  of  similarity  principles 
and  the  quasi-steady  assumption,  the  two-dimensional,  unsteady  gas-phase 
equations  have  been  reduced  to  a  one-dimensional  form.  Upon  integration  of 
that  ordinary  differential  equation,  one  function  results  that  can  he  input  in 
tabular  form  and  mav  be  considered  as  equivalent  to  an  algebraic  statement  and 


therefore  a  "zero-dinensional"  representation  of  the  gas-phase  has  Seen 
developed . 

The  liquid  hydrodynamic  field  is  represented  in  terms  of  the  classical 
Hill's  spherical  vortex  and  the  liquid  temperature  and  concentration  are 
assumed  to  be  uniform  (but  time-varying)  along  each  streamsurface .  This  high 
Reclet  number  approximation  allows  conduction  onlv  normal  to  the 
streamsurface.  The  characteristic  liquid  heating  length  and  heating  time  are 
thereby  significantly  reduced  from  the  stagnant  value. 

The  liquid  phase  heat  diffusion  equation  has  been  shown  to  be 
approximately 
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when  *  is  the  normalized  streamf unction  and  T  is  the  liquid  temperature.  Note 
that 


dr 


bo  R2 


dt 


and 


C  =  C(T)  = 
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The  boundary  conditions  on  (6)  are  given  as 


(i)  at  t  =  0,  T  =  0 


(ii)  at 
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In  the  above  eauations,  t  is  the  time  variable,  a  is  the  thermal  diffusivity, 
R  is  the  droplet  radius,  RQ  is  the  initial  radius,  o  ^  is  the  normalized  heat 
flux  at  the  surface,  and  b,  and  b  are  constants. 

On  account  of  C,  Equation  (6)  is  nonlinear.  However,  C  will  be 
considered  as  slowly  varying  compared  to  T  and  will  be  taken  as  piecewise 
constant  in  our  integration  calculation.  The  approximately  linear 
differential  operator  allows  us  to  develop  a  Green's  function  and  reduce  the 
solution  to  an  integral  form  whereby,  via  quadrature,  the  temperature  at  any 
point  can  be  determined  as  a  function  of  the  heat  flux  at  the  liquid 
surface.  In  particular,  an  integral  equation  can  result  that  relates  the 
surface  temperature  to  the  surface  flux.  The  Green's  function  or  kernel 
function  is  given  by 


K(t,*) 
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where  F  and  G  are  analytic  functions  of  4>.  The  eigenvalues  X  are  determined 

'  n 

by: 


G( >  ,C  ,1)  =  0  (in) 

n  o 

Now  the  temperature  anywhere  in  the  liquid  nay  be  determined  as  a  function  of 
the  surface  heat  flux 
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where 


f (  t)  =  yj  at  <fc  =  1 


In  particular,  at  the  surface,  <b  =  1,  we  find  that 
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From  the  gas-phase  analysis,  we  also  have  a  relationship  between  the 
surface  temperature  and  B'  where  B'  does  depend  upon  both  surface  temperature 
and  surface  heat  flux.  This  allows  for  complete  determination  of  the  surface 
temperature  and  heat  flux  upon  coupling  with  the  liauid-phase  solution.  The 
solutions  are  numerical  and  time-dependent. 

The  integral  equation  at  the  surface  can  be  replaced  bv  the  followine 
svstem  of  ordinary  differential  equations. 


T 
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I  g  (T) 
n=0  n 


First  we  define 
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Then  it  follows  that 


dg 

_ c 

dx 


=  a  f( x) 


and 


dg 


rr  ■  V„(T)  *  8„f(t>  !  "  > 


Note  that  by  definition 


g  (0)  =  0  for  all  n  >  0 
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(13) 
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The  independent  time  variable  in  Equations  (13)  and  (14)  is  a  Laeraneian 
time  so  that  an  advantage  of  these  equations  is  that  thev  can  be  solved  as  we 
follow  the  vaporizing  droplets  through  the  spray  combustor.  As  Figure  2 
shows,  less  computer  time  is  required  for  the  ordinarv  differential  equation 
svstem  solution  than  for  the  direct  numerical  solution  of  Equation  (6).  The 
improved  efficiency  measured  in  CP1T  time  occurs  at  anv  desired  accuracy. 

In  the  case  of  multicomponent  fuels,  liquid  mass  diffusion  is  of  critical 
importance.  The  more  volatile  components  will  vaporize  first  from  the  surface 
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leaving  the  less  volatile  components  as  the  doninant  species  in  the  vicinitv 


of  the  surface.  Liouid  mass  diffusion  brings  the  volatile  connonents  to  the 
surface  for  further  vaporization.  The  vaporization  rates  of  each  component 
will  change  with  time.  In  the  latter  part  of  the  droplet  lifetime,  tKe 
vaporization  rates  of  the  less  volatile  components  will  become  greater  than 
the  rates  for  the  other  components  whose  concentrations  has  decreased. 

It  is  necessary  to  solve  diffusion  equations  for  n-1  of  the  n  components 
in  the  liquid-phase.  Species  conservation  eauations  for  gas-phase  fuel  vapor 
components  must  also  be  solved.  The  same  type  of  algorithms  can  be  used  here 
as  were  used  in  the  solution  f  the  temperature  for  the  gas  and  the  linuid. 
While  interface  conditions  and  diffusivity  coefficients  differ,  the  same  type 
of  differential  operators  are  present.  Details  of  the  recommended  methodologv 
for  multicomponent  fuels  are  presented  in  References  11  and  12. 

The  results  of  Figure  3  show  surface  temperature  and  droplet  diameter 
versus  time.  All  quantities  are  normalized  and  nondimensional .  The  behavior 
is  clearly  seen  to  be  highly  transient  and  the  danger  of  oversimplified 
assumptions  about  surface  temperature  is  portrayed.  Figure  4  shows  the 
behavior  of  hexane  as  the  more  volatile  component  in  a  hexane-decane 
bicomponent  mixture.  A  general  decrease  in  surface  mass  fraction  and 
component  vaporization  rate  occurs  with  time  due  to  mass  diffusion  control.  A 
slight  increase  in  the  hexane  surface  mass  fraction  occurs  late  in  the 
lifetime  when  the  surface  temperature  becomes  large  enough  to  increase 
substantially  the  decane  vaporization  rate. 

Figure  5  shows  the  variation  of  the  coefficient  C  with  time.  Note  that 
during  the  earlier  part  of  the  droplet  lifetime,  the  assumption  of  slowlv- 
varving  C  appears  to  be  poor.  Nevertheless,  results  for  temperature, 
concentration,  vaporization  rate,  and  diameter,  such  as  those  shown  in  Figures 
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even  a  thousand  representative  droplets  must  he  considered 


Nevertheless . 


more  new  and  useful  information  does  result  from  more  detailed  calculations, 
especially  from  "exact"  finite-difference  computations  that  minimize  t^e 
number  of  approximations. 

The  simplified  method  uses  only  boundary  layer  theory  thereby  providing 
no  information  about  the  flow  aft  of  the  point  on  the  droplet  surface  where 
the  shear  stress  goes  to  zero.  The  drag  force  on  the  droplet  therefore  cannot 
be  calculated  from  that  method.  The  drag  force,  or  equivalent lv ,  the  drag 
coefficient  must  be  known  in  order  to  predict  droplet  trajectories  and  the 
local  vapor  mixture  ratio.  The  drag  coefficient  is  expected  to  he  verv 
different  from  that  of  an  isothermal  solid  sphere  due  to  a  combination  of 
effects.  Blowing  due  to  droplet  vaporization,  variable  properties  due  r  > 
temperature  and  molecular  weight  gradients,  apparent  mass  due  t a  flow 
acceleration  or  deceleration,  and  internal  liquid  circulation  will  eac^  i-pi.'t 
the  drag  coefficient.  Preliminary  results  from  a  number  of  investigators 
indicate  a  reduction  in  drag  coefficient  from  hard-sphere  results.  Mussel t 
and  Sherwood  numbers  also  experience  modification  from  hard-sphere  results. 

The  exact  solutions  are  also  useful  as  a  basis  for  comparison  for  the 
simplified  model.  Since  detailed  experimental  results  do  not  exist,  the  exact 
solutions  give  the  most  reliable  information  available  about  the  flow  field 
within  and  around  the  droplets.  The  validity  of  various  assumptions  used  in 
the  approximate  models  and  the  accuracy  of  those  models  can  be  evaluated 
through  the  use  of  the  exact  models,  since  the  exact  solutions  do  not  require 
any  limitation  on  Reynolds  number,  they  cover  the  zero  Reynolds  number 
(stagnant),  low  Reynolds  number,  and  high  Reynolds  number  (boundary  layer) 
domains  while  the  simplified  model  discussed  in  the  previous  section  is 
limited  only  to  the  high  Reynolds  nu'  her  domain.  Other  approximate  solutions 


exist  for  the  transient  heating  and  vaporization  in  the  stagnant  case. 


Dwver,  Sanders,  and  co-workers^  ,16,17,18,19)  ^ave  performed  the  ■one': 
work  on  this  subject.  In  a  series  of  papers^'  ^ ^  ^  ^  assuming  a  unif;>r~ 
densitv  gaseous  nedia,  thev  demonstrated  certain  interesting  results.  The- 
internal  circulation  resulted  in  vortical  streamsurfaces  that  did  rese-hle 
those  for  Hill's  spherical  vortex.  The  internal  circulation  did  affect  the 
transient  heating,  especially  at  higher  Reynolds  numbers,  and  transient 
heating  was  verv  important  to  the  droplet  vaporization  process. 

The  Reynolds  number  decreased  with  time  due  to  deceleration  of  the 
droplet.  For  solid  spheres,  the  drag  coefficient  increases  as  Reynolds  number 
decreases.  Their  results  with  vaporization  artificially  eliminated  do 
indicate  a  similar  increase  of  the  drag  coefficient  during  the  droplet 
lifetime.  However,  when  vaporization  occurred  in  accordance  with  the  natural 
laws,  the  drag  coefficient  decreased  with  time,  contrary  to  intuition.  It  is 
known  that  vaporization  reduces  friction  drag  due  to  the  boundary  layer 
blowing,  but  still  the  result  must  be  considered  as  remarkable  because  it 
constitutes  a  qualitative  difference  and  not  just  a  quantitative  difference 
when  compared  to  solid  sphere  data. 

The  variable  density  case  has  been  studied  by  Patnaik,  et  al^^)  and  the 
trend  for  the  drag  coefficient  is  the  same  generallv  but  even  more  dramatic 


than  in  the  constant  density  case.  Figure  7  shows  a  typical  result  whereby 
droplet  radius,  mass,  relative  velocity,  Reynolds  number  based  upon  relative 
velocity,  and  drag  coefficient  all  decrease  during  the  droplet  lifetime.  The 
rapid  decreases  of  velocity  and  Reynolds  numbers  are  especially  notable.  The 
vaporization  rate  first  increases  and  then  decreases  during  the  droplet 
lifetime.  The  Nusselt  number  decreases  during  the  droplet  lifetime  for  both 
the  constant  densitv  and  variable  densitv  cases. 


The  above  calculations  that  relate  to  decelerating  droplets  are  made  in  a 


frame  of  reference  fixed  to  the  droplet.  Since  such  a  reference  frame  is  no-'- 
Newtonian,  a  reversed  D'Alembert  force  must  be  applied  uniformly  throughout 
the  field.  That  is,  at  each  time-step,  there  is  a  change  in  relative  velocitv 
between  the  free-stream  and  the  droplet  that  occurs  because  of  the  ehanee  of 
velocitv  of  the  droplet  due  to  the  drag  force.  The  droplet  velocitv  is 
returned  to  zero  in  the  frame  of  reference  by  deceleration  of  the  flow 
field.  Mot  iust  the  free-stream  velocity  is  decreased  relative  to  the 
droplet,  but  is  decreased  relative  to  the  droplet,  but  the  whole  gas  flow 
field  is  reduced.  It  is  important  that  the  numerical  adjustment  of  the 
velocity  field  is  made  at  every  point  in  the  field  at  each  time  step  to 
properlv  represent  the  reversed  inertia  force.  It  is  clear  from  the  constant 
density  and  the  variable  density  calculations  that  the  drag  coefficient  will 
depend  upon  a  nondimensional  acceleration  or  deceleration  factor.  (The 
standard  parameter  is  an  acceleration  modulus.)  The  reversed  inertia  force 
will  make  the  drag  coefficient  different  from  its  steadv-state  value. 

A  general  observation  from  the  various  calculations  is  that  the  transport 
within  the  droplet  is  convective-dominated  in  the  earlv  portion  of  the  droplet 
lifetime.  In  the  later  part  of  the  lifetime,  the  droplet  relative  velocity 
and  the  strength  of  the  internal  vortex  have  been  reduced  substantially  so 
that  internal  conduction  has  increased  in  relative  importance.  This  indicates 
that  the  simplified  model  of  the  previous  section  can  be  appropriate  for  the 
earlv  lifetime  but  can  become  inaccurate  for  the  later  lifetime.  Clearly,  a 
more  general  simplified  model  is  reauired.  One  approach  is  to  use  the  vortex 
model  of  the  previous  section  for  a  portion  of  the  lifetime  and  then  to  switch 
to  a  spherically-symmetric  conduction  model  later  in  the  lifetime. 
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The  effect  of  burning  in  the  gas  film  and  wake  of  the  dronlet  was 

considered  bv  TVver  and  Sanders.^  ^  Thev  treated  an  uncounled  gas-phase 

problem  onlv  but  it  is  reasonable  to  believe  that  their  conclusions  will  not 

change  for  the  coupled  gas-liauid  problem.  Thev  found  a  thick  flame  raising 

2 

auestions  about  the  validity  of  a  d  law  for  stagnant  droplets.  For  Reynolds 
numbers  near  100,  thev  found  burning  only  in  the  wake  with  little  heat 
feedback  to  the  droplet. 

It  is  reasonable  to  believe  that  the  combined  effects  of  droplet  group 
behavior  and  high  Reynolds  number  will  limit  the  possibility  of  flames 
enveloping  individual  droplets  in  practical  combustors.  This  is  an  important 
issue  in  the  modelling  of  spray  combustors  since  we  must  know  whether  or  not 
to  include  the  heat  of  combustion  in  the  calculation  of  the  transfer  number 
F'.  Without  combustion  in  the  boundary  layer  or  near  wake  of  the  droplet,  the 

value  of  B'  would  be  lower. 

(18") 

Dwyer  and  Sanders  have  examined  the  effect  of  pressure  on  the  droplet 
vaporization,  on  the  drag  coefficient,  and  on  the  Nusselt  number.  They 
considered  three  cases  at  one,  five,  and  twenty-five  atmospheres, 
respectively.  Reynolds  numbers  and  Nusselt  numbers  decreased  with  time  in  all 
gases.  Using  a  normalized  time-scale  independent  of  the  pressure,  they  showed 
that  Reynolds  number  decreased  more  rapidlv  and  Nusselt  number  more  slowly 
with  time  as  pressure  increased.  Drag  coefficients  decreased  with  time  for 
one  and  five  atmospheres.  The  decrease  was  more  rapid  for  one  atmosphere.  At 
twentv-five  atmospheres,  the  drag  coefficient  increased  with  time. 

The  decrease  in  Reynolds  number  without  mass  transfer  has  been  shown  to 
result  in  an  increase  of  the  drag  coefficient  with  time.  Mass  transfer  tends 
to  reduce  the  drag  coefficient  and  mass  transfer  rates  increase  wit^  time 
during  the  droplet  lifetime.  Apparently  for  one  and  five  atmospheres,  the 


increase  In  mass  transfer  rates  more  than  compensates  the  decrease  In  Revnolds 
number.  Therefore,  the  draff  coefficient  decreases  with  time.  However  for  tKe 
largest  value  of  pressure  and  gas  densitv,  the  Reynolds  number  decrease  is  so 
rapid  that  even  with  increasing  mass  transfer  rates,  the  droplet  drag 
coefficient  increases  with  time. 

Tt  is  clear  that  an  extensive  parameter  survev  is  needed  to  determine  the 
drag  coefficient  as  a  function  of  the  instantaneous  Reynolds  number,  droplet 
acceleration  modulus,  and  the  instantaneous  transfer  number.  Such  numerical 
data  would  be  critical  input  to  a  spray  combustor  analysis.  The  preliminary 
evidence  cited  above  indicates  that  solid  sphere  data  overpredicts  the  drag 
coefficient  and  therebv  underpredicts  the  penetration  of  the  droplets  in  a 
sprav  combustor. 

Another  important  parameter,  especially  in  dense  sprays  involves  the 
effect  of  droplet  spacing  upon  the  droplet  vaporization  rate,  the  drag 
coefficient  and  the  Nusselt  number.  Preliminary  results^  '  indicate  that  the 
vaporization  rate,  Nusselt  number,  and  drag  coefficient  can  be  strongly 
affected  by  the  proximity  of  other  droplets.  The  flow  field  around  two 
droplets  in  an  axisymmetric  configuration  was  calculated.  Each  droplet  was 
artificially  kept  at  constant  Reynolds  number.  One  droplet  is  aft  of  the 
other  at  a  constant  spacing.  The  calculations  for  the  first  (upstream) 
droplet  are  made  neglecting  the  presence  of  the  second  (downstream)  droplet. 
The  wake  solution  of  the  first  droplet  then  provides  the  upstream  condition 
for  the  calculation  of  the  flow  around  the  second  droplet. 

Figure  8  shows  that  the  drag  coefficient  generally  decreases  with  time 
due  to  the  increasing  mass  transfer  effects.  The  drag  coefficient  is  higher 
at  the  lower  Reynolds  number  and  is  greater  for  the  first  droplet.  It  can  be 
expected  that  the  drag  coefficient  of  the  second  droplet  will  be  stronglv 


dependent  upon  the  distance  between  droplets.  Further  wor'<  on  the  ofrec:  of 
droplet  soacincr  is  required  in  the  future. 

Sone  interesting  calculations  are  presented  in  Reference  23.  The  effects 
o?  variable  properties  are  included  and  are  shown  to  be  of  some  importance, 
the  utilitv  of  the  analysis  is  limited  because  the  study  addresses  onlv  low 
values  of  the  transfer  number  and  the  reversed  inertia  force  is  not  oroperlv 
considered.  The  conclusions  that  drag  coefficient  increases  during  the 
droplet  lifetime  is  seriouslv  misleading. 

In  summary,  computational  fluid  dynamics  provides  a  useful  mechanism  for 
the  studv  of  fuel  droplet  vaporization  in  sprav  combustion.  Its  potential  has 
certainly  been  demonstrated  and  it  should  provide  valuable  inputs  on  drag 
coefficients,  Sherwood  numbers,  and  Nusselt  numbers  for  sprav  analyses.  The 
need  for  extensive  parameter  survevs  to  determine  the  dependence  of  drag 
coefficients  and  Nusselt  numbers  on  transfer  number,  Reynolds  number,  Prandtl 
number,  droplet  spacing,  and  droplet  acceleration  modulus  are  required. 

FLAMF  PROPAGATION  ANO  FLAMEHOLDING  STUDIES 


In  situations  where  droplet  heating  and  vaporization  are  rate-controlling 
and  ambient  temperatures  are  sufficiently  high  that  droplet  lifetime  is  not 
long  compared  to  droplet  heating  time,  the  details  of  the  transient  droplet 
thermal  behavior  are  very  Important.  It  has  been  shown  “  that  different 
vaporization  models  can  produce  s ieni f 1  cant  1 v  different  results  for  local 
gaseous  fuel-air  mixture  ratios.  Three  models  were  considered:  infinite- 
conductivity,  spherical lv-symmetri c  heat  conduction,  and  heat  conduction  with 
high  Peclet  number  internal  circulation  (vortex).  Differences  in  mixture 
ratio  were  as  high  as  inn?  especially  in  the  earlv  vaporization  period.  Those 
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differences  would  dramatically  affect  predictions  of  flame  stability 


a"  1 

ignition. 

On  account  of  the  significant  sensitivity,  it  is  necessary  to  employ  the 
most  accurate  model  in  a  given  combustion  situation.  For  zero  or  low  Fevnolds 
number  based  upon  relative  gas-droplet  velocity,  the  spherically-symmet ric 
droplet  heating  and  vaporization  model  is  recommended.  When  the  droplet 
Reynolds  number  is  large  compared  to  unity,  the  vortex  model  is  most 
appropriate.  Only  in  the  limit  where  droplet  lifetime  is  much  longer  than 
droplet  thermal  time  does  the  infinite  conductivity  model  give  accurate 
results.  Note  that  all  three  models  give  the  same  result  in  that  limit  and 
also  the  limit  is  not  interesting  in  combustion  applications. 

Eulerian-Lagrangian  methods^ ^  have  been  utilized  to  studv  flame 
propagation'  '  and  flame  holding'  '  in  various  configurations.  Transient 
droplet  heating  and  vaporization  were  considered  in  each  of  these  studies. 
Transient  droplet  behavior  was  also  taken  into  consideration  in  a  recent  paper 
on  ramjet  combustion  instability.'  ’ 

(25) 

A  one-dimensional,  unsteady  analysis  of  flame  initiation  and 
propagation  in  a  closed  volume  showed  manv  interesting  results.  The  droplet 
number  density  was  initially  uniform  but  expansion  of  the  burning  gases  and 
droplet  motion  due  to  drag  caused  subsenuent  stratification  of  the 
heterogeneous  mixture.  A  hot  wall  at  one  end  resulted  in  ignition  and  the 
propagation  of  a  flame  across  the  chamber.  The  fuel  droolets  do  not 
completely  vaporize  as  they  are  swept  by  the  flame.  This  is  especially  time 
for  larger  droplets  and  less  volatile  fuels.  The  droplets  vaporize  in  the  hot 
gas  behind  the  flame  front  and  fuel  vapor  diffuses  to  the  flame.  Figures  9, 
in,  and  11  show  that  the  advancing  flame  is  actually  a  diffusion  flane  with 


fuel  vapor  and  oxvgen  coming  from  opposite  sides  of  the  flame. 


Note  that  resolution  is  obtained  on  a  scale  smaller  than  the  distance 
between  droplets.  Figure  10  shows  local  peaks  of  fuel  vaoor  mass  fraction  at 
the  droplet  locations  which  are  source  points  for  the  fuel  vapor.  The  zone  of 
chemical  reaction  in  the  diffusion  flame  is  not  large  compared  to  droplet 
spacing  so  that  resolution  on  this  scale  is  absolutely  necessary. 

(  J 

A  two-dimensional,  steady  flame  propagation  analysis'  '  considers 
parallel  streams  of  droplets  flowing  adjacent  to  air  streams  that  initially 
are  cold.  The  domain  is  bounded  bv  two  walls,  a  hot  wall  and  a  cold  wall, 
with  heat  conducting  in  a  transverse  direction  across  the  flow.  Tn  the 
direction  of  flow,  vaporization  and  mixing  of  fuel  vapor  with  air  occurs 
first.  Then  the  combustible  mixture  ignites  and  as  indicated  by  the  schematic 
of  Figure  12,  a  premixed-like  flame  propagates  across  the  flow. 

All  of  the  droplets  are  not  vaporized  in  this  first  flame-front;  rather, 
vaporization  continues  beyond  the  front.  All  of  the  oxygen  is  not  consumed  in 
the  first  flame  front,  either.  Diffusion  flames  are  established  around  the 
droplet  streams  in  the  hot,  oxygen-rich  gases  downstream  of  the  premixed-like 
flame  as  shown  in  Figure  12.  Fuel  vapor  mass  fraction  contours,  shown  in 
Figure  13,  support  the  schematic  description  of  Figure  12. 

The  results  of  References  25  and  26  together  show  that  the  simple  concept 
of  a  flame  sweeping  a  heterogeneous  mixture  and  resulting  in  conversion  of 
liauid  fuel  to  gaseous  products  in  a  short  distance  is  verv  seriously 
flawed.  Except  for  verv  small  droplets  and  verv  volatile  liquids,  the  flame 


asvroototieal  ly  in  time  via  finite-difference  techniques.  polvdi sper<~  an! 

nono-’isperse  soravs  have  been  considered  in  the  calculations.  Transient 

droplet  heating  and  vaporization  are  considered  following  the  methodology 

outlined  earlier  in  this  paper.  4  Lagrangian  calculation  is  made  for  droplet 

properties  and  an  F.ulerian  calculation  is  employed  for  the  gas 
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properties.  '  A  standard  k,  e  turbulence  model  is  utilized. 

Figure  14  shows  droplet  trajectories  and  size  for  a  oolydisperse  sprav 
calculation  while  Figure  15  indicates  gas  velocity  magnitude  and  direction  at 
the  grid  points  of  a  nonuniform  finite-difference  mesh.  It  is  seen  that  a 
double-vortex  recirculation  zone  occurs  behind  the  center-body;  due  to  the 
momentum  of  the  liquid  sprav  originating  from  the  center-body  at  the  axis  of 
symmetry ,  the  inner  vortex  results.  The  smaller  droplets  in  the  spray  tend  to 
become  bound  in  the  recirculation  zone  while  the  larger  droplets  penetrate 
bevond  the  recirculation  zone.  Note  that  solid  sphere  drag  data  is  enploved, 
thereby  overpredicting  the  drag  coefficient  and  underpredicting  the  droplet 
penetration.  In  these  calculations,  the  recirculation  zone  was  predicted  to 
be  verv  fuel  rich  and  turbulent  transport  of  fuel  vapors  from  the  fuel  vapor 
from  the  recirculation  zone  to  the  incoming  air  was  also  predicted  to  be  very 
important.  Improved  drag  coefficient  calculations  should  decrease  the 
importance  of  turbulent  transport  of  fuel  vapors  and  increase  the  relative 
importance  of  droplet  penetration. 

concluding  REMARKS 


An  integrated  approach  to  spray  combustion  modelling  requires  several 
tvpes  of  analyses  as  key  elements.  Analysis  of  individual  droplet  behavior 
with  the  intention  of  developing  an  algorithm  or  a  sub-grid  model  for 
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vaporisation  is  required. 


An  analysis  to  determine  drag  coefficients 


for  u  s  l- 


in  droplet  traiectorv  analvsi?  is  mandatorv.  Finallv,  a  finite-difference 
analvsis  of  the  sprav  combustor,  emploving  inputs  from  the  droplet  analyses, 
is  renuired. 

Numericallv-ef f icient  and  phvsicallv-accurate  algorithms  for  transient 
droplet  heating  and  vaporization  have  been  developed.  Effects  of  internal 
liquid  circulation  have  been  carefully  considered  where  appropriate.  These 
models  have  successfully  been  applied  to  combustor  calculations. 

Much  work  on  the  exact  solution  of  the  flowfields  within  and  surrounding 
droplets  has  been  performed.  Drag  coefficients  for  vaporization  an i 
decelerating  droplets  have  been  shown  to  be  very  different  from  those  for 
isothermal,  constant  velocity,  solid  spheres.  A  more  extensive  parameter 
studv  is  still  required  to  provide  useful  input  to  spray  combustor  analyses. 
It  should  show  the  dependence  of  drag  coefficient  and  Nusselt  number  upon  mass 
transfer  number,  acceleration  moduli,  and  inter-droplet  spacing. 

Spray  combustion  analyses  have  been  performed  and  results  presented  for 
one-dimensional  and  two-dimensional,  steady  and  unsteady  flows  with  and 
without  recirculation  and  with  monodisperse  and  polvdisperse  spravs. 
Lagrangian  droplet  calculations  and  Eulerian  gas-phase  calculations  have  been 
made.  The  need  for  improved  droplet  drag  coefficient  data  is  a  maior 
limitation.  Another  limitation  is  the  need  for  a  dronlet  vaporization  model 
that  applies  over  a  wide  range  of  Revnolds  numbers. 
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Figure  1.  Vaporizing  Convective  Droplet 
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Figure  2.  Maximum  Percentage  Error  in  Mass  Vaporization 
Rate  Versus  CPU  Time  Comparisons 


25 


TRANSFORMED  TIME 


Figure  5.  Negative  Coefficient  C  Versus  Time  Comparisons 
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Figure  6.  Surface  Temperature  Versus  Time  Comparisons 
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Figure  7.  Overall  Droplet  Behavior 
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Figure  8.  Drag  Coefficient 
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Figure  9.  Gas  Temperature  VS  Axial  Position 
at  Different  Times 
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Figure  10.  Fuel  Vapor  Mass  Fraction  VS 

Axial  Position  at  Different  Times 


AXIAL  POSITION  (CM) 


Figure  11.  Oxygen  Mass  Fraction  VS  Axial 
Position  at  Different  Times 


Figure  12.  Vaporization,  Ignition,  and  Figure  13.  Fuel  Mass  Fraction 

Combustion  of  Two  Parallel  Contours 

Fuel  Droplet  Streams 
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Figure  14.  Droplet  Trajectories  in  Polydisperse 
Spray  with  Recirculating  Flow 
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Figure  15.  Gas  Velocity  in  Recirculating  Flow 
with  Polydisperse  Spray 
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